On Collective Effects in Cavity Quantum Electrodynamics 
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I. INTRODUCTION 



The idealized system of a two-level atom interact- 
ing with a second quantized single-mode electromagnetic 
field, confined in a cavity, plays an important role in the 
study of various fundamental aspects of quantum me- 
chanics. The micromaser is a remarkable experimental 
realization of such a simple but fundamental system (for 
reviews and references see e.g. 0). It is therefore also 
an example of one of the rare systems in Nature which 
exhibit a rich structure of physics that can be investi- 
gated experimentally and which, at the same time, can 
be studied by exact theoretical methods. In the optical 
regime a microlaser has also been realized experimentally 
U Recently this has also been achieved for a one-atom 
system Q. 

Many features of the micromaser system can be re- 
garded to be of general interest. Various aspects of 
stochastic resonance has e.g. been explored in this 
system The micromaser also illustrates a feature of 
non-linear dynamical systems: turning on randomness 
may led to an increase of the signal to noise ratio 0, @ . 
It can, furthermore, be argued that the micromaser sys- 
tem is a simple illustration of the conjectured topological 
origin of second-order phase transitions [g, |jj . Trapping 
states H, U have been generated in the stationary state 
of the micromaser system and therefore the generation 
of states with no classical analogue in such a system is 
feasible [lOl [Tl| . A basic ingredient of the micromaser 
is the description of the dynamics in terms of the Dickc 
model fl2| in the so called rotating wave approximation, 
i.e. the of the Jaynes-Cummings (JC) model [lijj. Early 
experimental studies involves a confirmation of the JC- 
model predicted atom revivals ,14]. Recently one has 
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also explicitely demonstrated field mode quantization in 
a cavity . Entanglement of mesoscopic states of the 
electromagnetic field and an atom in a cavity has also 
been demonstrated ^(| in accordance with theoretical 
considerations (see e.g. Ref. 0])- Even though our 
analysis will focus on the dynamics of the micromaser 
we observe that the JC model with d amp ing effects in- 
cluded has been realized in ion traps |18j and in super- 
conducting systems |l9j . The coupling of electromagnetic 
modes to the latter artificial two-level systems has been 
demonstrated in the laboratory [2fl l2l| . We also notice 
that it has experimentally been verified that phonons can 
be confined in semiconductor planar cavities p3 |. It is 
therefore not unlikely that the present analysis may find 
applications in systems similar to the micromaser system 
but realized in a completely different physical framework. 

The paper is organized as follows. In Section [H] we 
outline the dynamics of a typical experimental setup of 
the two-level system interacting with a single-mode of 
the radiation field. Long-time correlations are discussed 
in Section IIIII and corrections due to detection efficien- 
cies are discussed in Section IIVI Collective effects due 
to the finite probability of having two atoms at a time 
in the cavity are analyzed in Section [V] together with 
possible effects on the detection of trapping states. In 
Section I VII we summarize our work and indicate effects 
on the phase transitions of the micromaser system in the 
so called large- TV limit. 



II. MICROMASER DYNAMICS 



In our analysis we consider the following typical real- 
ization of the micromaser. The pump atoms which enter 
the cavity are at resonance with the radiation field of the 
cavity, and are also assumed to be prepared in the ex- 
cited state. The injection intervals between the incoming 
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atoms are assumed to be Poisson-distributed. In terms 
of the dimensionless atomic flux parameter N = R/"/, 
where R is the rate injected atoms and 7 is the damping 
rate of the cavity, the stationary photon number proba- 
bility distribution is then described by a diagonal density 
matrix with diagonal elements which are well known ||| 
and are given by 



Pn = Po 



n 



n b m + Nq m 
(1 + rib) m 



(1) 



Here q m = q(x) = sin 2 (9y r x), with x — m/N, and where 
we have defined the natural dimensionless pump param- 
eter 9 = gry/N in terms of the atomic transit time r. 
Furthermore, g is the single photon Rabi frequency at 
zero detuning of the JC-model ^| • The overall constant 
po is determined by Y^=oP™ = 1- 

The theory as developed in Refs.|3.l23j suggests the ex- 
istence of various phase transitions in the large-TV limit 
as the parameter 9 is increased. A natural order pa- 
rameter is then the average photon "density" (x) , where 
( ) denotes an average with respect to the distribution 
Eq. (JJJ. An exact large- N limit treatment of the mi- 
cromaser phases structure and the corresponding crit- 
ical fluctuations in terms of a conventional correlation 
length was presented in Refs. 24]. Spontaneous jumps in 
(n) /N and large correlation lengths close to micromaser 
phase transitions have actually been observed experimen- 
tally nm. 

Several new intriguing physical properties 
of the micromaser system are unfolded when the theo- 
retical analysis is extended to a more general setup of 
the parameters available in the micromaser system than 
those considered here 261. 



the JC-model QjJ |24| . The lowest eigenvalue Ao = of L 
then determines the stationary equilibrium solution p = 
as given by Eq. (jTJ. T s (t) = Tr[M(+)p] = u T M(s)p 
then is the probability that an atom is found in the state 
s = ±, where +(— ) denotes the excited(ground) state, af- 
ter it leaves the microcavity. The vector u T is the trans- 
pose of the vector u with all entries equal to 1. Here we 
have used the fact that the vector u T simply represents 
the trace operation and p is the diagonal density matrix 
of the cavity field. When the injection intervals between 
the incoming atoms are Poisson-distributed, as we have 
assumed is the case, the joint probability, V(s\, s 2 ,t), of 
observing two atoms in the states si and S2 with a large 
time delay t between them can be written in the form 



V(s u s 2 ,t) = Tr[5( S2 )e-^*5(s 1 )p] 
= u T M{s 2 )e^ Lt S{ Sl )p 



(4) 



The time delay t corresponds to a large number, k ~ 
Rt, of unobserved atoms between the two detections. In 
Eq. @ we make use of the propagation matrix 



S(s) = (l + L c /N)- 1 M(s) 



(•5) 



where S = S(+)+S(—). S is a so called stochastic matrix 
(see e.g. Refs.|M EU) with left and right (p^) 

eigenvectors corresponding to the eigenvalue K n , where 
n = 0,1,2,..., such that 1 = kq > K\ > ... > 0. The 
stationary distribution as given by Eq. corresponds to 



the eigenvalue kq = 1, i.e. vffl = u and 
spectral decomposition of S , where k = 0, 1, 



„(o) 



p. The 
i.e. 



(6) 



III. LONG-TIME CORRELATION EFFECTS 



Let us now consider long-time correlations in the large- 
N limit as was first introduced in Refs.p^|. These cor- 
relations are most conveniently expressed in terms of the 
continuous-time formulation of the micromaser system 
|27| . The vector p formed by the diagonal density matrix 
elements of the photon field then obeys the differential 
equation 



dp 



(2) 



where L = Lc — N(M—1). Here Lc describes the damp- 
ing of the cavity , i.e. 

(Lc)nm = (rib + 1) [ nS n>m - (n + l)S n +i, m ] 

+ n b [(n + l)5„, m - n<5„, m+ i ] , (3) 



and M = M(+) + M(-), where M(+), 



q n +i)S n , m and M(- 



is useful in many of the numerical calculations presented 
below. Here the normalization is such that u^ T p^ 1 — 
S nm . The joint probability V(si, S2,t) is symmetric, i.e. 
V(s\,S2,t) = V(s2, si,t), and is properly normalized, i.e. 
Ssis 2 ^ , ( Sl ' s 2j^) — 1 since Sp = p and u T p = 1. In the 
original discrete formulation of the micromaser system Q 
the corresponding joint probability Vk(si,S2) of observ- 
ing two atoms in the states si and S2, with k unobserved 
atoms between, can be written in the form 

Vk(si,s 2 ) = Tr[S(s 2 )S k S( Sl )p] = u T M(s 2 )S k S( Sl )p . 

This joint probability is also properly normalized, i.e. 
Ssis 2 'Pkisi, S2) = 1. For a sufficiently large k and with 
k ~ Rt it now follows that 

V k (s 1 ,S2)=u T S(s 2 )e-7 Lt S(s 1 )p = V(s 1 ,S2,t) (8) 

in the large- N limit [24[ and where L is as in Eq. @ ■ A 
formal way to see the validity of Eq. (JSJ is to notice that 
(1 + Lc/N)- 1 = 1 - L c /N + 0(1/N 2 ). We can then 
write 



qnfin,m+i have their origin in 



5 = 1 



N 



(L c - N(M - 1)) 



(9) 
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apart from terms of order Lq{M — 1)/N and higher or- 
ders in l/N. Below we will find expansions of the form 
Eq.© very useful in the analysis of detection efficien- 
cies. We observe that the joint probabilities Vk(s\,S2) 
are symmetric in s\ and S2 |24j as in the continuous- 
time formulation, i.e. the order in which the atoms are 
measured is irrelevant. Other definitions for such joint 
probabilities have appeared in the literature. In the def- 
inition for the two-time coincidence probability used in 
e.g. Refs. |30H3l| . S(s 1 ) in Eqs. (QJ and J7J is replaced by 
M(si). With the Poisson-statistics assumption behind 
the derivation of these equations this would, however, 
seem like an unnatural thing to do. Since the correlation 
length to be defined and used below is only sensitive to 
the eigenvalues of the operator L, our discussions below 
will not, in the end, be effected by such a modification, at 
least not in the large- N limit. With detector efficiencies 
■q + and r\- for detecting atoms in the excited or ground 
state respectively, the sequence probabilities defined and 
used in Ref.|3^| are identical to our joint probabilities 
Eq.Q in the discrete formulation of the micromaser if 
i] + = 7]- = 100%. For k — the joint probability 
'Pk{si,S2) reduces to T , SiS2 (t), i.e. the probability that 
the next atom is in the state S2 = ± if the previous 
atom has been found in the state si. V Sl (r) exhibits the 
experimentally observed revivals in microcavity systems 
and in e.g. ion-traps T~ > s 1 s 2 ( T ) exhibits in 

addition so called pre-revivals |24j. 

A properly normalized correlation function 7c (t) can 
now be defined and expressed in different but equivalent 



manners, i.e. 



7c(t) 



1 - (s) 2 



V(-,-,t)-V(-) 
V(+)V(-) 



V(+,+,t)-V(+) 2 
V{+)V{-) 

r(+)v(-)-r(+,-,t) 



V(+)V(-) 



(10) 



where (ss) t = J2 Sl , S2 sis 2 V(s 1 , s 2 , t) and (s) = 
^2 s sV(s). This correlation function satisfies —1 < 
lc{t) < 1- At large times t — ► 00, we then define the 
atomic beam correlation length £c by [24| 



7c (t) = ice * /€c , 



(11) 



which then is determined by the next-to-lowest eigen- 
value Ai of L, i.e. 7^(7 = 1/Ai. In the discrete formula- 
tion of the micromaser system we make the replacement 
V(si,S2,t) — > Vk{si,S2) in Eq.lJTDJ and, for a sufficiently 
large k, we then define, in a similar manner, the correla- 
tion length £d by the expression 



lD(k) = 7ce 



(12) 



In the large-TV limit one can show that £c and £d con- 
verge to the same limit and we therefore write £ = — 
£d for sufficiently large TV [2J|. For photons a similar 
analysis leads us to a correlation length £ 7 . It follows 
that the correlation lengths so defined are identical in 
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FIG. 1: The convergence of the correlation length 
1/.R£d(/c) = — log(|7r>(fc)|)/fc as a function of the number 
k ~ Rt of unobserved atoms leaving the cavity for 9 — 0.5, 
iV = 10 and rib ~ 0.15 in the case of the discrete formulation 
of the micromaser system with an analytical fit using only the 
next to the leading eigenvalue ki of the stochastic matrix S 
as defined in the main text. The corresponding correlation 
length is R£ w 14. The upper curve corresponds to the def- 
inition of joint probabilities of the present paper. The lower 
curve corresponds to a redefinition S(si) — > M(si) in the 
joint probabilities Eq.JTJ. The horizontal lines corresponds 
-log(«„) for n = l, 7. 



the large- N limit, i.e. = ^ 0- In Fig. [2 we illus- 
trated the convergence of 1/R£n(k) = — log(|7c(A;)|)/fc 
to its asymptotic value 1/RS, = — log(/«i) for a typical 
experimental setup of the micromaser system. We ob- 
serve that spectral resolution Eq.JfJJ and the definition 
Ea.lfTOjl modified to yield 7n(fc) leads to 



lD{k) 



where we have defined 



, exp(-felog(l/K„)) 



(13) 



■S1S2 



u T M(+)p { s n) u { s n)T S(+)p 
V+V- 



(14) 



The analytical fits in Fig^ ^/R^d — ci/fc — log(/€i), are 
based on using only the n = 1 term in Ea. (|13(l . We have 
found that this approximation works for all 6 despite the 
fact that N is small in this case. The numerical value of 
Ci can easily be computed numerically given the micro- 
maser parameters. The upper curve in Fig^corrcsponds 
the definition Eq.Q of joint probabilities and leads to 
Ci w 3.40. The lower curve in Fig. ^ corresponds to 
a redefinition S(s±) — > M(s\) in the joint probabilities 
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Eq.Q and leads to c\ w 4.37. In Fig[3]we compare the 
correlation length 7£ using the discrete and continuous- 
time formulation of the micromaser system for a typical 
setup of parameters. We observe the rapid convergence 
of the two formalisms already for a small value of N. In 
the sequel we will therefore make use of the fact that 
7£ = 1/Ai = l/iVlog(l//ci) if N is large enough. 
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FIG. 2: The correlation length 7^ as a function of 8 for N — 10 
and rib = 0.15 in the case of the continuous-time (solid line) 
and the discrete (dashed line) formulation of the micromaser 
system as described in the main text. 



IV. DETECTION EFFICIENCIES 



It has been emphasized in the literature that in the 
analysis of the micromaser system, when detection ef- 
ficiencies has to be taken into account, one must dis- 
tinguish between the occurrence of a certain detection 
event and the absence of such a detection event. This 
is so since the detection of an atom leaving the cavity 
necessarily gives rise to an instantaneous and non-local 
state reduction of the micromaser cavity radiation field. 
In Refs. [32] a non-linear equation of motion was derived 
taking such effects into account. In studying the statis- 
tics of sequence events it was, however, shown in Ref . |3dj| 
that one may make use of a linear equation of motion 
for non-normalized states. Recently it has been shown 
in detail how these different formalisms actually leads to 
the same physical results [33| . Below we follow the anal- 
ysis of Refs. [3(1 l33l] since we find it simpler to implement 
numerically. 

Our calculation of joint or sequence probabilities when 
detection efficiencies must be taken into account follows 



the scheme outlined in Ref. [33j . There is simple prescrip- 
tion how to modify the analysis in Section liTTl to the sit- 
uation with, in general, different detection efficiencies 77+ 
and r\-. The basic matrices M(±) are naturally modified 
to M{s) = T] s M(s) with s = ± and we therfore also mod- 
ify M = M{+) + M(-) to M = f]+M{+) + n-M(-). It 
is also convenient to define the matrix M_ = M — M. It 
then follows that the joint probability as given by Eq.Q) 
is modified according to 



Pk(si, s 2 ) -> Pk{si, a 2 ) = u T M(s 2 )S k S{ Sl )p/M 



where we have defined the matrix 



S(s) = (1 + L c /N + M_)- x M(s) 



(15) 



(16) 



Here S = S + + S- and TV = i] + V+ + r/-V- is a 
normalization factor. We also make the modification 
V{s) — > V(s) = r) s V(s)/J\f. It is of importance to observe 
that the stationary micromaser photon number distribu- 
tion Eq. (JJJ still can be obtained from the stationary con- 
dition Sp = p. The discussions in Section IlIII in the case 
of the discrete formulation of the micromaser system can 
now simply be carried through by changing the relevant 
probabilities by the modified probabilities defined above 
and one, e.g., shows that "Pk( s li s 2) — ~Pk(s9., s-\), which 
was shown explicitely for k = in Ref. [33] • We therefore 
define the correlation function 



lD(k) = 



?(+)?(-) -7M+,-) 
P(+)P(-) 



(17) 



expressed in one of many equivalent manners as in 
Eq. (|17|) . For N large enough we can then write £ = £d if 
i] + = Tj- = 1. Using a 1/N expansion similar to Eq.@ 
we now find that 



NM_) + {M-l) = S 



(18) 



apart from higher order terms in (Lq — NMJ)/N and 
higher orders in 1/N. For sufficiently large N we would 
then e.g. conclude that 



Pk(si,s 2 ) = u T M{s 2 )e^ Lt S{ Sl )p/N 



(19) 



The correlation length £ = £d in Ea. l|17|) would then 
again be determined by the next-to-lowest eigenvalue 
of L and, as a result, £ would actually not depend on 
the detection efficiencies at all and therefore equal to £. 
This was the claim as announced in Refs. 26]. This con- 
clusion is basically correct apart from a trivial correc- 
tion which, unfortunately, was not taken into account. 
The physical explanation of this correction is simply 
the fact that the time interval t in Eq. l|19|l should be 
such that Rt is average number of atoms which the ob- 
server claims passes through the micromaser system, i.e. 
Rt ~ k(T] + V(+) + r)-V{-)) and not Rt ~ k. This 
means that the correlation length £ is " renormalized" to 
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£ = (rj + 'P(+) + r)-V(—))£. A more formal proof of this 
assertion can be given as follows. Let us first consider 
the case of equal detection efficiencies, i.e. r\ = rj + = r)-. 
Since the correlation length £ will be determined by 
the next-to-leading eigenvalue R\ of the operator S, i.e. 
7£ = l/iVlog(l/rei) for large N, we study the eigenvalue 
problem SpD = RpD- As in Ref. |24j it is convenient to 
rewrite such an eigenvalue problem in an equivalent form, 
i.e. 

(h c - N(l + ? - n)(M - 1))) p D = r)N(l - l)p D ■ 

( 20 ) 

This equation can now be compared to the eigenvalue 
problem in the case of the continuous-time formulation 
of the micromaser system with 77^ = 77 = 1, i.e. 



(L c - N{M - l)) Pc = X(N)p c 



(21) 



A(jV(l+?- 77))= 77iV(~-l) 



(22) 



1 

— = 1 



Ai_ 

r]N 



1 



1 



(23) 




where we have made the N dependence explicit in the 
eigenvalue X(N). By comparing Eas. (|20() and (|21|l we 10 
conclude that 



N = 10 
I n b = 0.054 

ii 

- rj + V++rj-V-„ 



Since \(N) remains finite in the large-^V limit, we can 
write A (iV(l + -q/R — r/)) — » X(N) for sufficiently large 
N and we conclude that 



We therefore find that 7^ = 777^ apart from 1/N correc- 
tions. A consequence of this analysis is that r)N(l/R — 
l)Mpr> = T]N(1/k— 1)Pd in the large- N limit as can 
be seen again by comparing Eqs.J^OJ) and (|21|l . This 
can be understood as follows. For large TV the com- 
ponents eigenvector pu will be distributed around some 
large component. We can the replace the matrices M(s) 
for s = ± by the diagonal matrices 'P(s), i.e. their mean- 
field values, in the expression r]N(l/R — l)MpD above. 
Extending the analysis above for a equal detection ef- 
ficiencies to the case with different ones, we then find 
that 7<£ = (rj + V(+) + rj^V(~ ))7£ apart from 1/N correc- 
tions. Numerically, it turns out that the convergence of 
the correlation length as a function of N is very rapid. In 
Figs. |31 we study, as an example, the correlation length 
7^ = j£d for a moderate value of N = 10 and other 
parameters adapted to experimental data on trapping 
effects in the micromaser system [HI . One verifies that 
7^(6* = 0) = (./Vlog(l + l/iV?7 + ))~ 1 in excellent agreement 
with numerical simulations in general. We find it re- 
markable that the scaling law for detection efficiencies of 
the ratio of correlation lengths £/£ = r/ + V(+) + Tj-V(— ) 
works so well for such a small value of N as used in e.g. 
Fig s El Experimentally one could therefore measure the 
appropriate sequence probabilities, evaluate the correla- 
tion length and then renormalize the corresponding data 
with an easily calculable factor. The resulting correlation 
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FIG. 3: In the upper graph the correlation length 7^ is given 
as a function of 6 for N = 10 and rif, = 0.054 for various val- 
ues of detection efficiencies. With equal detection efficiencies 
■q = 77+ = rj- (lower curve) the correlation length is given by 
7/7^ with £ = = r i- = !)■ I n the lower graph we con- 

sider the same parameters but a " renormalized" correlation 
yt/{v+'P+ + V-'P-)- 



length so obtained is then the one that was predicted in 
Refs.li. 



COLLECTIVE EFFECTS 



Collective effects are now due to the fact that during 
a time interval t, such that < t < r, two atoms has 
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jointly interacted with the same cavity radiation field. 
An ideal and very special situation corresponds t = r 
which has been discussed in great detail in Ref.[24| and 
also elsewhere [35| . The general situation is more tedious 
but straightforward to analyze and has been discussed in 
great detail in Ref. 36] in terms of an expansion in the 
parameter e = Rt. This parameter is supposed to be 
small, i.e. e <C 1, in order to be close to the one-atom 
maser situation (see e.g. Appendix A in the second ref- 
erence of Ref. [24[ ) • We have reconsidered the analysis of 
Ref. [3(| . The generator L in the master equation Eq. J2J) 
is, up to first order terms in e, modified according to 

fr= - j(L c -(l-2e)N(M-l))p , 

+ jN— f dtu 2 {t)p , 



and, finally, 
d n {t) = 4 



n + 1 n + 2 2 
2n + 3 q 

n + 2 



2„ - :',2n + 3 q ^ m - qn+l{T - t)] 

q n +3/2(2t)q n+1 (T - t) . (30) 



2(2n + 3) 

Our final expression for L tot is, in fact, in agreement with 



(24) 

where the two-atom generator u 2 (t) describes two atoms 
that have jointly interacted with the cavity during the 
time interval t such that < t < r and given explicitely 
in Ref. [3(| . One finds that the generator L in the master 
equation Eq.(0| is replaced by L to t = L + L co i, where 
L co i describes the two-atom collective effects with matrix 
elements given by 

(L col ) nm = Ne[(w n (T) + v n (r))5 

n,m 




(25) 



Here we have defined 
( \ 1 r 

V n [T) = 



dt ( [q n+ i(T -t) + &„(*)] [1 - q n +l(T - t) 

- qn+2(r - t)] - g„+i(r)[l - q n +i{r) - gn+2(r)] 

- [c„,(t) + d n (t)][q n+ i(r - 1) - q n+2 (T - t)]J , 

(26) 



and 



w n {r) = i J dt(^C n (t) + q n+2 {T - t)q n+1 (T - t) 
+ b n (t)q n+2 (T -t) + [c n (t) + d n (t)]q n+2 (r - t) ) 




(27) 



In these definitions we make use of the functions q n (t) 
q n = sin 2 (gty/n), 



n + 1 

K{t) = 2n q n+ 3/ 2 (2t)[l - g ra+ i(r - t)] 
- ^q n +3/2(2t)q n +i(T - i) , 



(28) 



and 



n + 1 

C ™W = olo^ i o\ 1™+3/2(2t)[l - gn+l(r - t)\ 



2(2n + 3) 
+ ll+s^Qn+iir - t) , 



(29) 



FIG. 4: In the upper graph the order parameter (a;) = (n/N) 
is given as a function of 9 for N — 50 and nt — for various 
values of e = Rt. It is only for e = that trapping-state 
peaks are visible. In the lower graph the correlation length 
7£ is given for the same set of parameters. Even though the 
trapping-state peaks vanishes as e increases, the correlation 
length can still be large for e 7^ 0. The peaks at 6 ~ 1 
correspond to the first maser transition. 

the result of Ref.j3|J using a slightly different notation. 
It is rather straightforward to implement the matrix el- 
ements of L to t hi a numerical routine in order to find 
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the new stationary distribution p corresponding to the 
eigenvalue A = of L tot to be used in evaluating vari- 
ous expectation values. The correlation length is again 
determined by the next-to- lowest eigenvalue Ai of L to t- 
In the numerical work it turns out to be sufficient to use 
200x200 matrices for both L and Ltot in order to obtain 
the accuracy of the graphs as presented in the present pa- 
per. In Fig0]we show the results of a numerical evalua- 
tion of the order parameter (a;) and the correlation length 
7^ for N = 50 in a vacuum configuration (nj, = 0). 



atom inversion J(r) is given by 



TV = 10 
n b = 0.001 



e = 0.00 
e = 0.05 
e = 0.10 
e = 0.50 




J(r) = 



- v-V(-) 

TV 

+ n>b V++V- 



N 



(31) 



where is M is the normalization factor rj + 'P(+)+r]^'P(— ). 
Eq. (|31|) gives the general relation between the atomic in- 
version 7(r) and the order parameter (x). We observe 
that if 77 = r/ + = r/- then J(t) is independent of r\. It 
is clear form our Fig|B] that trapping effects are much 
more visible in the correlation length than in the order 
parameter (x) . Its is also clear that collective effects are 
small with the set of micromaser parameters chosen. In- 
stead detection efficiency will be a major correction to 
the theoretical values. In Section HVI we have seen how 
the correction length is to be corrected for due to de- 
tection efficiencies. In view of these results it therefore 
appears that trapping-states can be more clearly revealed 
experimentally in terms of the correlation length rather 
than the atomic inversions or, equivalently, the order pa- 
rameter (x). 



FIG. 5: The correlation length 7^ as a function of 9 for various 
values of e. The vertical lines indicate the trappings values of 
9 = tvVN(1/VS, 1/V2, 1, 2/V3, V2, VS, 2). 



Our results basically agree with the corresponding nu- 
merical results of Ref.|36( even though we appearently 
have a higher numerical precision. As argued in Ref . |36| . 
effects of trapping states at 8 = kn^N/n, for k, n — 
1,2,..., vanishes as e 7^ 0. This is, however, not so for 
the correlation length 7^ as indicated in the lower figure 
of Fig In fact, if we consider even a lower value of 
N = 10 and n b = 0.001 as in FigEI trapping-state ef- 
fects are clearly seen in the correlation length. In Fig [5] 
the vertical lines corresponds to the clearly visible trap- 
ping states with e — 0. In the experimental study of 
trapping-states in Ref.[lfJ the parameters are, however, 
varied in such a way that the parameter e is not constant 
but is given by e = Rt = Q^fN^/g as in Fig|SJ where 
we have chosen the cavity temperature to be T = 0.3K, 
corresponding to n/, = 0.054. The average photon life- 
time in the cavity corresponding to Figgis 0.1 s, i.e. 
7 = 10 8 . The vertical lines in Fig|5] correspond to the 
trapping-states as considered in Ref. [Kg and studied in 
terms of the atomic inversion. With detection efficien- 
cies taken into account, the appropriate definition of the 




FIG. 6: The order parameter (x) and the correlation length 
7^/5 as a function of 9 with no collective effects taken into 
account (e = 0) or e = Rt = OyN^f/g. The Rabi fre- 
quency is as in the experimental study of trapping-states in 
Ref-Efll, i.e. g = 39 kHz and 7 = 10 s" 1 . The vertical 
lines indicate the experimentally observed trappings values of 
9 = k^/N(1/V3,1/V2,2/V6,2/VE,1,2/^,V2) as reported 
in Ref. PH. 
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VI. FINAL REMARKS 



In conclusion, we have studied how to take detection 
efficiencies into account when comparing experimental 
data for the correlation length to the theoretical predic- 
tions. We have found a remarkable and simple scaling 
relation that connects observational data and the theo- 
retical prediction. We have also studied two-atom col- 
lective effects in terms of the natural parameter e = Rt. 
Even though cumbersome, the calculations are in prin- 
ciple straightforward. Trapping effects appear to sup- 
pressed with increasing values of e. We have, however, 
seen that in a realistic experimental situation, as dis- 
cussed in e.g. Ref . [Iol| . collective effects are, neverthe- 
less, small and detection efficiencies are more important 
to take into account. As we have shown, detection effi- 
ciencies can, however, be taken into account in a straight- 
forward manner. 

With increasing values of N, signals due to the micro- 
maser phase transitions become more pronounced. Gen- 
eral methods, which are exact in the large-iV limit, for 
computing characteristic features of these phase transi- 
tion have been presented in Rcfs. 24, 26]. Detection effi- 
ciencies will only imply a calculable rescaling as discussed 
in Section llVl 



i* a* a* a* 

^01 w 12 w 23 




5 10 15 20 



FIG. 7: The correlation length 7^ as a function of 9 for various 
values of e. The vertical lines indicate the large N values for 
the maser phase transitions at 6 = 9g — 1, #01 ~ 6.6610, # 12 ~ 
12.035, #23 ~ 17.413. With (x) as an order parameter, the 
transition at 6 = 8q is second order while the other are first 
order transitions. 

Due to the random arrival statistics of the pump atoms 



there is, however, in an actual experimental situation a 
finite probability of more than one pump atom in the 
cavity 34]. Since the average number of pump atoms 
inside the cavity is e = tR = ^^/Nd/g, this parameter, 
as we have seen above, naturally parameterize the prob- 
ability of collective pump atom effects. If e is assumed to 
be constant and small, i.e. e -C 1, we see that the dimen- 
sionless pump parameter N is bounded by y/~N -C g/j&- 
Arbitrarily large values of N can then only be reached 
by making g/j arbitrarily large which, of course, is dif- 
ficult to achieve in a real experimental realization of the 
micromaser. Alternative e = 6\fNg/^ should be small 
as 9 varies. As we have seen in the present paper, if e is 
sufficiently small, corrections to the observables as dis- 
cussed in the present paper can be calculated in a rather 
straightforward manner. 

At finite N and including collective pump atom ef- 
fects, signals of the large N phase transitions in the 
order parameter (x) are still clearly exhibited, at least 
in the case when the pump atoms are prepared in the 
excited state and at resonance with cavity radiation 
field. The critical point 0q of the first second-order 
maser transition remains the same with (x) used as a 
natural order parameter. As seen from Eq. l|24[l when 
taking collective effects into account for a sufficiently 
small e, one-atom effects are modified by a renormal- 
ization N — > exp(— 2e)N, where exp(— 2e) w 1 — 2e is 
the probability for one-atom events in the micromaser 
cavity. We therefore expect that the critical parame- 
ters 0* kk+1 cx "JN of first-order transitions are changed to 
^fc/c+i = exp (— e)6Zu ji ■ By making use of the general re- 
sults of Refs.[24l26j concerning the N dependence of the 
peak values of the correlation length these are changed 
accordingly, and we find that \og{"l £,) C rit is changed to 
Iog(70crtt = exp(-e)log( 7 £(0 = 9* kk+1 )). In FigQwe il- 
lustrate these effects in the case N = 100 and rib = 0.15. 
For small values of e we find that the scaling behavior 
given above for 9 kk+1 and log(7£) cr i t are indeed compat- 
ible with the exact results. 

As in Ref. [3^] we also find that our perturbative meth- 
ods gives meaningful results even when the expansion 
parameter e is large. In the present paper we have consid- 
ered a limited range of micromaser parameters. If atoms 
and the radiation field are not at resonance one can e.g. 
extend the results of the present paper by making use of 
the methods of Ref. [2^] and the procedures outlined in 
the present. 
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